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Abstract. The motion of a collisionless plasma - a high-tcmpcrature, low- 
density, ionized gas - is described by the Vlasov-Maxwell system. In the pres- 
ence of large velocities, relativistic corrections are meaningful, and when sym- 
metry of the particle densities is assumed this formally becomes the relativistic 
Vlasov-Poisson system. These equations are considered in one space dimension 
and two momentum dimensions in both the monocharged (i.e., single species 
of ion) and neutral cases. The behavior of solutions to these systems is studied 
for large times, yielding estimates on the growth of particle momenta and a 
lower bound, uniform-in-time, on norms of the charge density. We also present 
similar results in the same dimensional settings for the classical Vlasov-Poisson 
system, which excludes relativistic effects. 

1. Introduction. The fundamental equations which describe the time evolution of 
a collisionless plasma in the presence of large velocities are given by the relativistic 
Vlasov-Maxwell system (RVM). The simplest form of these equations that retains 
electromagnetic effects can be obtained by posing the problem in one space dimen- 
sion and two momentum dimensions, the so-called "one and one-half" dimensional 
system with a single species of ion: 

dtf + vAf + (E 1 + v 2 B)d v J + {E 2 - viB)d v J = 

p(t,x)= / f(t,x,v)dv, j(t,x)= / vf dv 



Ei{t, x ) = \(^j v) d v - J p(t, v) d v 

d t B + d x E 2 = 0, d t E 2 + X B = -j 2 



(RVM) 
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where 



is the relativistic velocity. If one further assumes that the field components E2 and 
B are initially zero and the number density / is initially even in v%, then these 
qualities persist for all t > 0, and the system reduces to the "one and one-half" 
dimensional relativistic Vlasov-Poisson system: 

d t f + vid x f + E 1 d Vl f = 



p(t,x) 



E x {t,x) 



f(t,x,v) dv 
P(t,y)dy- / p{t,y)dy 



(RVP) 



Here, it should be noted that / = f(t,x,vi,v 2 ) depends on two components of 
momentum. Finally, if one wishes to consider a classical system without relativistic 
effects, then the velocity is represented by v <= K 2 and not v. In this way, one then 
obtains the "one and one-half" dimensional classical Vlasov-Poisson system: 

d t f + v 1 d x f + E 1 d v J = 



Ei(t,x) 



p(t,x) = j f(t,x,v) dv 

2 ( / P(t,y)dy- j p{t,y)dy 



(VP) 



Though the difference between (RVP) and (VP) seems quite subtle (merely a "hat" 
on v is included or excluded), it is well-known that the behavior of these two sys- 
tems can drastically differ. This will be further highlighted by Theorem 1.1 which 
displays a significant difference in the long-time dynamics of the charge density 
when comparing the classical and relativistic systems. In each of the systems above 
t > represents time, x G K is position, v S K 2 represents momentum, E and B 
are electric and magnetic fields, / is the number density in phase space of particles, 
and all physical constants, including the particle mass and speed of light, have been 
normalized. 

In addition to the monocharged problems (RVP) and (VP), we will consider the 
analogous neutral problems that include an arbitrary number N species of ion, given 
by: 

dtf a + vxd x f a + e a E x d v J a = 



p(t,x) = J ^2e a f a (t,x,v) dv 



Ei(t,x) 



and 



Eh.(t,x) 



P{t,y)dy- / p(t,y)dy 

I J X 

dtfa + vid x f a + e a Eid Vl f a = 
p(t,x) = j ' s ^2,e a f a {t,x,v) dv 

•* a 

1 



(RVPN) 



p(i, y)dy 



p{t, y)dy 



(VPN) 
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respectively, where we have a Vlasov equation for the number density f a of each 
species indexed by a = 1, ...,N, and e a is the charge of the a species. Again, the 
f a depend upon two components of momentum v± and V2, and physical constants, 
such as particle mass, have been normalized. In the cases of (RVPN) and (VPN), 
we will assume the condition of neutrality 



for all t > 0. Another main focus of the paper will be to demonstrate the differences 
in qualitative behavior which arise when one compares monocharged and neutral 
plasmas, such as (RVP) with (RVPN) or (VP) when compared with (VPN). This 
is displayed by the differences in large time behavior of particle momenta for the 
systems (Theorems 1.3, 1.4, 1.5, and 1.6). 

Over the past twenty years, considerable progress has been made regarding the 
well-posedness of classical solutions to the Cauchy problem for Vlasov-Maxwell and 
Vlasov-Poisson set in a variety of dimensions (see [6] ) . Though the issue of smooth 
global existence for arbitrary data remains an open question in three dimensions 
for both relativistic problems, it has been proven for many lower-dimensional ana- 
logues, such as (RVM) (see [7]), as well as, for the classical Vlasov-Poisson system 
posed in three dimensions ([13], [15], and [18]). More specifically, it is well known 
that solutions of (RVP), (RVPN), (VP), and (VPN) remain smooth for all t > 
with f(t, •, •), or f a (t, •, •) for multi-species problems, compactly supported for all 
t > 0, assuming that the data possess the same properties. While these exis- 
tence theorems have contributed greatly to the mathematical understanding of the 
equations, very few results concerning time asymptotic behavior have appeared in 
the literature. Due to the so-called "dilation identity" , some time decay is known 
for Vlasov-Poisson in the classical, three-dimensional case ([10], [12], [14]). Addi- 
tionally, there are time decay results for the monocharged plasma (VP) when / is 
independent of V2 ([1], [2], [17]). The only known results regarding long-time behav- 
ior for a relativistic plasma were obtained in [11] and, more recently, [8]. These deal 
with the three-dimensional, monocharged problem with spherical symmetry and the 
one-dimensional, neutral problem with two species, respectively. References [3], [4], 
and [5] are also mentioned since they deal with time-dependent rescalings and time 
decay for other kinetic equations. We cite [16] as a general reference regarding the 
study of the Vlasov-Poisson system. 

Due to the large assortment of phenomena that may be exhibited by plasma of 
different species, one expects the behavior to be quite complicated and depend on 
many factors, including the size of the total charge, the sign of the net charge, and 
the variety of ionic species that are involved. For simplicity, we focus on cases in 
which the plasma is monocharged (i.e., composed of a single species of ion) or is 
composed of an arbitrary number of species and satisfies the condition of neutrality 
(i.e., possesses zero net charge). Hence, the methods used in the previously men- 
tioned articles do not apply. The present work, then, seeks to determine information 
regarding the large time behavior of solutions to (RVP), (VP), (RVPN), and (VPN) 




(N) 



a 



which, by conservation of charge, guarantees 




p(0, x) dx = 
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under these conditions. We remark that all of the results which follow will continue 
to hold in the strictly one-dimensional case, in which the number density / (or f a 
if multiple species are considered) is independent of v 2 - However, we keep the v 2 
dependence throughout because solutions to these equations also satisfy (RVM), 
thereby providing some information regarding time asymptotics in this case. In 
addition, a previous result [8] concerning a neutral plasma in one dimension did not 
readily generalize to the case in which v 2 is included. Thus, we believe it necessary 
to present the more general results for the "one and one-half" dimensional prob- 
lems, rather than asking the reader to believe that one-dimensional results could 
be easily extended to this case. Of course, we are also interested in the long-time 
behavior of (RVM), and this problem is explored in the recent paper [9]. 

Throughout the paper, we make the assumption that the initial data /(0, x, v) = 
fo(x,v), or in the case of multiple species f a (0,x,v) = fo a (x,v), are C 1 and com- 
pactly supported. We first consider the relativistic, monocharged case (RVP) and 
state the main result that, unlike solutions in the classical case, those which satisfy 
(RVP) do not give rise to a charge density that decays in time. 

Theorem 1.1. Consider (RVP) and assume fo G C^(M 3 ) is not identically equal 
to zero, then there exists C > such that 



for all t > 0,p G [1, oo]. 

Since [8] showed some time decay for solutions to (RVPN), in the case N = 2, 
e a = ( — 1) Q+1 , and f a independent of v 2 , this theorem displays the distinct dif- 
ference in behavior between monocharged and neutral plasmas, even in a one- 
dimensional setting. Hence, the dispersive effects which occur in the neutral case 
seem to be a stronger or more influential phenomena than the repulsion of the 
monocharged case. Additionally, while [1] showed that the charge density decays in 
sup-norm like t" 1 for the classical, monocharged system (VP) with / independent 
of V2 (and the argument can be extended to include solutions of (VP)), we find 
here that the charge density does not decay in any LP norm for the relativistic, 
monocharged system (RVP). Thus, we have differing asymptotic behavior depend- 
ing upon the inclusion or exclusion of relativistic velocity corrections. This can 
be contrasted with Horst's discovery [11] that solutions to both the classical and 
relativistic 3D Vlasov-Poisson systems satisfy the same time decay estimates under 
the assumption of spherical symmetry 

For these systems, the time asymptotic behavior of solutions will depend strongly 
on characteristics. Since we are first interested in (RVP), define the associated 
characteristics X(s,t,x,v) and V(s,t,x,v) by 



IIp(*)II p >c 




^=E 1 (s,X(s)) 
as 

V 2 (s,t,x,v) = v 2 



} 



(1) 



X(t, t, x, v) = x 
Vi(t,t,x,v) = vi 



LONG-TIME BEHAVIOR OF PLASMA IN 1.5D 



5 



Then, in addition to lower bounds on the charge density, we can also determine a 
uniform upper bound under a very general assumption. 

Theorem 1.2. Consider (RVP) and assume there is Fq 6 L 1 (M 2 ) such that 
fo(x, v) < F (v) Vx 6l,t)£l 2 . Then, 

dVi 

1 < —±(s,t,x,v) forO<s<t 
ovi 

and 

p(t, x) < J F (v) dv forf>0,j;e M. 

For each of these systems it is well-known that ||/o(£)||i is constant due to charge 
conservation. Hence, using Theorems 1.1 and 1.2 and interpolation with ||/o(f)||i, we 
may conclude in the case of (RVP) that ||p(i)|| P is O(l) for any p e [1, oo]. Finally, 
we can also sharply determine the asymptotic behavior of particle momenta for 
large time. 

Theorem 1.3. Let f(t,x,v) be a solution of (RVP) and define 

Qi(t) = sup{|wi| : there are x, V2 G K s.t. f(t, x, v) ^ 0}. 
Then, there are T, Ci, C% > such that for t > T, we have 

Cit < Q 1 (t) < C 2 t. 

To further display differences in behavior between monocharged and neutral 
plasma, wc consider the multi-species problem (RVPN). Since attractive forces 
among particles are now introduced, one intuitively expects that the particle mo- 
menta are asymptotically slowed in comparison to those of (RVP). This is accurate 
and demonstrated by the following theorem. 

Theorem 1.4. Let f a (t, x, v) satisfy (RVPN) for a = 1, ...,7V and define 
Qi(t) = sup{|«i| : there are x,i>2 G R s.t. f a (t,x,v) ^ 0}. 

Then, there is C > such that for any t > we have 

Qi(t) < cVTTt. 

Next, we turn our attention to the classical, monocharged problem (VP). As 
opposed to (RVP), time decay of the charge density is known from [1], so one may 
expect faster growth of particle momenta, as well. However, since the electric field 
terms are the same, the momenta are shown here to grow at exactly the same rate 
as in the relativistic case. 

Theorem 1.5. Let f(t,x,v) be a solution of (VP) and define 

Qi(t) = sup{|wi| : there are x, v-i € K s.t. f(t, x, v) ^ 0}. 
Then, there are T,C\,Ci > such that for t > T, we have 

C l t<Q l {t)<C 2 t. 

Finally, to contrast the result for (VP), we may show that particles travel at 
slower speeds in the neutral case (VPN). In fact, using very different techniques 
than the proof of Theorem 1.4, we obtain the same bound on growth as for the 
momenta in (RVPN). 
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Theorem 1.6. Let f a (t,x,v) satisfy (VPN) for a = 1, ...,N and define 
Qi(t) — sup{|ui| : there are x, v 2 G R s.t. f a {t, x, v) 7^ 0}. 

ot 

Then, there is C > such that for any t > we have 

Qi(t) < cVTTt. 

The next section is devoted to proving Theorems 1.1, 1.2, and 1.3. The proofs of 
Theorems 1.3 and 1.5 will be combined since they use similar techniques. Section 
3 then contains the proof of Theorem 1.4, while Section 4 contains the proof of 
Theorem 1.6. Throughout the paper, "C" will denote a generic constant which 
may change from line to line and depend upon initial data, but not on t, x, or 
v. However, constants with subscripts (e.g., "Co") will denote the same numerical 
value. 



2. Behavior of solutions to (RVP). Throughout this section, we will make 
great use of the monocharge assumption, which implies p > and thus d x E > 0. 
In addition, we will frequently use the quantity 

M := J J f {x,v)dvdx > 0. 

In order to prove Theorem 1.1, we first need a few lemmas. The following result 
shows that characteristics in (1) display an order-preserving property due to the 
monocharge assumption. 

Lemma 2.1. Let x < x* , v\ < v\, \v%\ > |«a |> an d < t\ < t^. Then, 

X(t 2 ,t u x,v) < X(t 2 ,h,x*,v*) 

and 

Vi(t 2 ,ti,x,v) < V\(t 2 ,ti,x* ,v*). 

Also, x < x* and t\ < t 2 implies X(t 2 , t\, x, v) < X(t 2 , t±, X* , v*), and similarly for 
characteristic momenta Vi if V\ < v\ and t\ <t 2 . 

Proof. Suppose x < x* , ti < t 2 , and define 

t = sup{t e [ti,t 2 ] : X(s,ti,X,v) < X(s,ti,x*,v*) Vs G [h,t]}. 

Then, for s G [ti,r], since p > 0, Ei(t, ■) is increasing and 

^-[V 1 {s,t 1 ,x\v*)-V 1 {s,t 1 ,x,v)] = £?!(«, ^(a.ti, !*,«•)) 
ds 

- Ex(s,X(s,h,x,v)) 
> 0. 



Thus, 



Vi(s, ti, x* , v*) — Vi(s, ti,x,v) > v\ — Vi > 
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and since ii\ = . Vl is increasing as a function of vi, we find 
V 1+ M 2 



X(s,ti,x*,v*) 



Vi(s,ti,x*,v*) 



> 



> 



yi + ^(Mi,a;*,v*) 2 + (^) 2 

V(s,ti,x,v) 

^l + Vi(s,ti,x,v) 2 + (v* 2 ) 2 

V(s,ti,x,v) 

y/l + Vi{s,t u x,v) 2 + V 2 . 
X(s, h,x, v). 



Hence, 



X(s, ti, x*, v*) — X(s, ti, X, v) > x* — x > 
and it follows that r = t 2 . A similar argument can be used to show the conclusion 
for characteristic momenta when Vi < v\. The cases x* = x or v* = v follow by 
continuity with respect to initial conditions. □ 

We will be concerned with extremal values of position and momentum on the 
support of / and will make great use of Lemma 2.1, so define the quantities 

W 2 = sup{|w 2 | : Bx, Vi e R s.t. f (x, v) ^ 0}, 



s(t) = G M. 2 : 3v 2 G R s.t. f(t,x,v) + 0}, 

and 

S{t) = s(t) x [-W 2 ,W 2 }. 
Next, we define the largest position and v\ momentum in this set, so let 
Pi(t) = max{vi : 3x,f 2 e R s.t. (x,v) e S(t)}, 
R(t) = max{x : 3v G R 2 s.t. (x, v) G S(t)}, 

and 

pi(t) = max{«i : 3v 2 G R s.t. (R(t),v) G 5(t)}. 

From (1) we know that V%{s,t, x, vi,V2) = v 2 for any choice of s,t,x, and v\. 
Thus, the ^-support of f(t,x,v) is constant in time. Hence, let p 2 = so that 
\ v 2\ > \p 2 \ f° r every v 2 with (x, u) G S(t). For the sake of notation, we will write 
P 2 {t) = P 2 = p 2 (t) = p 2 for every t > 0. Then, from the momentum functions, 
we define the vectors P(t) = (Pi(t), P 2 (t)) and p(t) = (pi(t),p 2 (t)), and finally the 
position 

r(t) = max{x : (x,P(t)) G S(t)}. 
Notice that these definitions concern the largest positions and V\ momenta, in ad- 
dition to the smallest values of \v 2 \ on the support of /, and thus, will allow us to 
utilize Lemma 2.1. For intuitive purposes, these quantities are illustrated roughly 
in Figure 1 below, which suppresses the dependence on v 2 . 

Since the field Ei only attains values in [— \M, |M], we have the following lemma. 
Lemma 2.2. For t 2 >t x > 0, 

Pi(ta)-Pi(ti)>^Af(t 2 -ti) (2) 



R{t 2 )-R(h) = f pi(s) ds. 

Jt! 



8 



ROBERT GLASSEY, STEPHEN PANKAVICH AND JACK SCHAEFFER 



Pl(t) 





^Wt),Pi(t)) 








I (R(t),pi(t)) 


\ s(t) 






X 


R(t) 



Figure 1. Maximal positions and momenta: An example of the 
support s(t). 

Proof. Let (X(t),V(t)) be characteristics of / as in (1). Partition [ti,t2] into n 
subintervals t\ = so < s\ < ■ ■ ■ < s n = ti and let 

Sfc+i - Sfc = As — 

n 

for all k = 0, n — 1. Then, for any k 

X(sk+i, Sk+i, R(sk+i),p(sk+i)) = R(sk+i) > X(sk+i, Sk, R(sk),p(sk)) (3) 

and hence the function £(t) := X(t, Sk+i, R{s k+ i),p(s k+1 )) - X(t, s fe , R(s k ),p(s k )) 
satisfies f(sfe + i) > 0. Similarly, 

X(sk,s k+1 ,R(s k+1 ),p(s k+ i)) < R{s k ) = X(s k ,s k ,R(s k ),p(s k )) 

and thus £,(s k ) < 0. So, by the Intermediate Value Theorem there exists r € 
[fife, Sfe+i] such that £(t) = and 

X(t, s k+1 , R(s k+1 ),p(s k +i)) = X(t, s fe , i?(sfe),p(s fe )). 

If we suppose 

Vi(r, s k ,R(s k ),p(s k )) > Vi{t, s k+1 , R(s k+1 ),p(s k+1 )) 
then by Lemma 2.1 and since 

\V 2 (t, s k ,R(s k ),p(s k ))\ = \p2(s k )\ = \p 2 (s k+1 )\ = \V 2 (t, sfc+i, R(s k+1 ),p(s k+1 ))\ 
we would have 

X(s,s k ,R(s k ),p(s k )) > X(s, s k +i,R(sk+i), p(s k+1 )) 
for all s > r, which contradicts (3). Hence, we find 

Vi(r, s k ,R(sk),p(s k )) < Vi(T,Sk+i,R(s k+ i),p(s k +i)). 
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By Taylor's theorem, 

Vi(T,* fc ,.R(* fc ),p(* fc )) > pi(s fe ) + J B 1 (s fc ,J?(s fe ))(r-s fc )-C(As) : 
= Pi(.s fe ) + iA/(r-.s fc )-C(A S ) 2 



and 



Vi(T,8k+l,R(Sk+l),p(Sk+l)) < Pl{Sk+l) + Ei(s k +l,R{sk+l)){T - Sfc+i) 

+ C(As) 2 

= Pi(«fc+i) + o M ( T _ + c '( As ) 2 - 

Hence, 

< Vi(t, s fc+ i,i?(s fc +i),p(sfc+i)) - Vi(r, s k , R(s k ),p(s k )) 
< Pi(s k+1 ) + ^M{T-s k+1 ) + C(As) 2 

Pl (s k ) + lM(T-s k )-C(Asf 

= pi(s k+1 ) -pi(sfe) - -M(s fc+1 - s fe ) + C(As) 2 

and we find 

Tl-l 



Pi(*2) -Pi(ii) = X! bi( s fe+i) --Pi( s fc)] 

1 



fc=0 

Tl-l 

> E 

fc=0 



^M( Sfe+1 - s fe ) - C(As) 2 



This holds for any n £ N, and thus it follows that 

Pi(t 2 ) - Pi(h) > -M(t 2 -h). 

In addition, p\ and p\ are both increasing and hence integrable. Next, 
R(s k +i) - R(s k ) „ . . 

Pl(Sfc) 

Sfe+i — Sfe 

X(sfe+i,Sfe,i?(Sfe),p(sfe)) -R(s k ) 

> Pl(SfcJ 

Sfe+l - Sfc 

_ X(s k +i,s k , R(s k ),p(s k )) - X(gfc, Sfe, R(s k ),p(s k )) 

Sfe+l - s fc 

- X(s k , s k , R(s k ),p(s k )) 
> -CAs 
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and similarly 

R(s k+1 ) - R(s k ) 



■Pi(sfc+i) 



Sfc+l ^ Sfc 

R(sk+i) - X(sk, sk+i, R(s k +i),p(sk+i)) 



< 



Pl(Sfc+l) 



Sfe+l - s k 

_ X(sk+i,Sk+i,R(sk+i),p(sk+i)) - X(s kl Sk+i 1 R(sk+i),p(sk+i)) 

Sfe+l - s k 

- X(sk+i,s k+ i,R(sk+i),p(sk+i)) 
< CAs. 

Multiplying these inequalities by Sk+i — Sk and summing over k = 0, n — 1, we 
find 

n-l 

fl(t a ) - JZ(*i) - ^pi(s fc )(sfc+i - s fe ) > -C(t 2 - h)As 

fe=0 

and 

n-l 

i?(i 2 ) - JJ(t!) - ^p 1 (s fc+1 )(s fe+1 - 8 k ) < C{t 2 - h)As. 

k=0 

Noting the missing t 2 term in the sum of the first inequality and the missing t\ term 
in the second, we find 



k=0 



< CAs + \pi(ti)\As + \pi(t 2 )\As. 



Finally, this holds for any n € N, and it follows that 



R(t 2 ) - R(ti 



pi(s) ds. 



Remark 1. The result (2) yields 

Pi(t) > pi (t) >pi(0) + ]-Mt. 
We may proceed in a similar manner to determine the behavior of r and P\ . 
Lemma 2.3. For t 2 > t\ > 0, we have 



□ 



(4) 



r(t 2 )-r(h) > / P 1 (s) ds, 

Jt! 

Ei{t 2l r{t 2 )) >Bi(ti,r(ti)), 

and 

Pi(t 2 ) - Pi(h) = ( 2 ^i(s,r(s)) ds. 

Jt! 

Proof. Proceeding as in the previous lemma, 

Pi{sk+i) > V 1 (s k+ i,Sk,r(sk),P{sk)) 

and 

Pi{sk) > Vi(sk,Sk + i,r(s k +i),P(sk+i))- 



(5) 
(6) 

(7) 
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Thus, there is r € [sfe,Sfc+i] such that 

Vi(t, s k+ i,r(s k+1 ),P(s k +i)) = Vi{t, Sk,r(s k ),P(s k )). 
Using Lemma 2.1 again and noticing 

\V 2 (r,s k ,r(sk),P(s k ))\ = |P 2 (s fe )l = l^(*fc+i)l = l^(r, a k+1 ,r(s k+1 ), P(s k+ i))\, 
we find 

X(t, s k+1 ,r(s k+1 ),P( Sk+1 )) > X(t, s k ,r(s k ),P(s k )) 
follows immediately. Now, 

X(t, s k+1 ,r(s k+1 ),P(s k+1 )) < r(s k+1 ) + A(sk+i)(r - s k+1 ) + C(As) 2 

and 

X(t, s k , r(s k ), P(s k )) > r(s fc ) + Pi(s*)(t - s k ) - C(As) 2 , 

so 

< r(s k+1 ) - r(s fe ) - P 1 {s k+1 ){s k+1 - r) - Pi(s,.)(r - s k ) + C(As) 2 
and hence 

n-i 

r(i a ) - r(* x ) > ^ [£i(*k+i)(*k+i - r) + Pi(s fe )(r - s fc )] - Cn(As) 2 . 

k=0 

Using the Lipschitz continuity of P\ , we find 



fe=0 



r(t 2 ) - r(ti) > 51 A(*fc)(*fc+i - - Cn(As) 2 
and (5) follows. Next, let 

f(s)=r(h)+ I H(r)dr. 



From the definition of -P(t) and as in the proof of Lemma 2.1, we have for (x, v) € 
S(t), 

Since ^(^2) > ^(^2), and i?i (t, •) is increasing, we find 



Ex(t 2 ,r(t 2 )) - Etih^ih)) > E x (t 2 ,f{t 2 )) - Exit^fih)) 



t'2 



dt 

(t,f(t)) 



f(t,f(t),v) (PlW-f'l 



dvdt 



> 

and (6) holds. Next, 

Pi(s fe+ i) - Pi(s fe ) 



Ei(s k ,r(s k )) 
Sk+i - s k 

> Vi(s k +i, Sk,r(s k ), Pjsk)) - Vi{sk,s k ,r(s k ),P(sk)) 

Sfe+l - s k 

~ Vi(s k , Sfe,r(sfc), P(sfe)) 

> -CAs 
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and similarly 

Pi(sk+i) - Pi{sk) „ t , u , r<A 

- Ei{s k+ i,r(s k+ i)) < CAs. 



Sfe+l - s k 



Hence, 



and 



n-l 



Pi(t 2 ) - Pi(h) - J2 Ei(*k,r(a k ))(a k+1 - s k ) > -Cn(As) 2 



k=0 
n-l 



Pi(t 2 ) - Pi(h) - J2 Ei(sk+i,r(s k+1 )){s k+1 - s fe ) < Cn(As) 2 



k=0 



Using (6) and the boundedness of E\ yields 



n-l 



Pi(t 2 ) - Pi(h) - J2 Ei(s k+1 ,r(s k+1 ))(s k+1 - s fc ) 



< Cn(As) 2 + CAs 



for every n £ N and (7) follows. □ 

Corollary 1. Ei(s,r(s)) — > — M as s — > oo. 

Proof. From the previous results (4) and (7), we may deduce 

Px(0)+ / £?i(«,r(«))dfl = i\(t)>pi(t)>pi(0) + -Mt 



o 



and hence for all t > 0, 



£i(s,r(s)) ds> -Mt-C. 



Additionally, we know that £a(s,r(s)) is increasing and E x {s,r{s)) < \M for all 
s > 0. Hence, Ei(s,r(s)) must converge as s — ► oo and its limit can be no greater 
and no less than \M . □ 

Finally, we may prove the non-decay result. In the proof of Theorem 1.1 the 
generic constant "C" may depend upon the fixed variable t > 0. 



1 f x 1 f°° f x 1 

E\ = - I pdy - - pdy= pdy - -M 



Proof of Theorem 1.1. Note that we may rewrite E\ as 
\j Pdy ~\ 

^ -J — oo ^ - . 

so that 

= -M + E x (t,x). 

By the corollary, we may choose t large enough so that 

£i(i,r(t)) > 0AM 

and hence 

r(t) r rPi(t) 

/ / f(t, x, v)dv\dv 2 dx > 0.9M. 

-oo J J — oo 

By the Intermediate Value Theorem, there is e > such that 

r r{t)-e r r Pi(t)-e 

f(t, x, v)dv 1 dv 2 dx = 0.8M. 
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Additionally, notice that 

r (t) r r Pi(t) 



f(t,x,v)dvidv2<ix > 

lr{t)-eJ JP 1 {t)-c 

since (r(i),P(t)) <G S(t). If (X(t),V(t)) is any characteristic then 

d r X(s) ^ fViis) 

f(s, x, v)dvidv2<ix 



ds 



— oo 



> o 

Similarly, we find 

We take- 
and 

Then, from above 
and 



d_ 

ds 



f(s,x, Vi(s),v 2 ) (E 1 (s,X(s)) - Ei(s,x)) dv 2 dx 
f(s,X(s),v) (V^s) - urj dvidvi 



OO p pOO 

I I f(s 1 x 1 v)dv\dv2dx > 0. 

X{s) J JVi{ S ) 



OO /> /"OO 



X( S )=X( S ,i,r(i)-e,P 1 (i)-e,P 2 ), 
Vi (s) = Vi (s, f , r(t) - e, P x (f ) - e, P 2 ), 

V 2 (a) =P 2 (s) = P 2 . 

r/ .V!( S ) 

/(s, x, v)dv 1 dv 2 dx > 0.8M (8) 

f(s,x,v)dvidv2dx>C>0 (9) 

for s > t. From (8) and the definition of Pi it follows that 

Pi(s,X(s)) > 0.3M 

and 

Vi(s) > 0.3M(s-t). 
dh 

Since /i(i>i,w 2 ) = ui satisfies — — = (1 + I ^ 1 2 ) 3 / 2 ( 1 + v \) > and is thus increasing 

OVi 

in V\ , we find for s > t 

Us)> , ■ (10) 

y/l + [0.3Af(s - <)] 2 + Pf 

Put A = 0.3M(s — i). From (10) and the compact w 2 -support of / we find 



1-ViM < 



VI + ^4 2 + P 2 



2 



1+P 9 2 



y/i + A 2 + P'i \Jl + A 2 + P 2 2 + A 
C 



~ l + A 2 
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and thus 



1-Vi(s) < 



C 



(11) 



1 + [0.3M(s - t)] 2 

for s > t. Because |Vi(s)| < 1, we see from the compact x-support of / that 
/ (t, x, v) = for x > R(0) + 1. Thus, p(t, x) = for x > R(0) + t. Finally, (9) yields 



C < I p(s, x) dx 

X(s) 
R(0)+s 

p{s, x) dx 

X(s) 

< Wp^W^RW + s-Xis)) 1 -*. 

However, using (11) 

R(0) + s-X(s) = R(0)+t-X(t) + (s-t)-(X(s)-X(t)) 

= R(0)+t-X(t) + J (l-Vi(r))dT 

C 

t 1 + [0.3Af(r - t)]- 



< C + 

< C 

where C may depend upon the fixed time t but not on s > t. So, we find 

0<\\ P (s)\\ P 
and the proof of the theorem is complete. 



□ 



Proof of Theorem 1.2. We will sometimes use an abbreviated notation for charac- 
teristics, namely 

X(s) = X(s,t,x,v) 
and similarly for V(s). Let t > and define 

OX 

t* = inf{r E [0, t) : — (s, t, x,v)<0 Vs E (r, t)}. 
ovi 



Note that (1) implies 

d fax 

ds \dvi 

and 



= (i+M 2 r 3/2 (i+, 2 2 )>o 



dX 



= 



dX 



so that t* is well-defined. Next, for every s G [r*,i] we know — — (s) < so 



and thus 

Using this lower bound, 
d fdX\ 



d_ ( dV\\ _ dX < 



1 + vk dVi 



= 1. 



> 



l + «2 



ds V^i/ (i + |y| 2 ) 3 / 2 avr - (i + |u| 2 ) 3 / 



>o, 
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and hence 

— (r*,t,x,v) < 0. 

dVi 

From this, it follows that t* = and hence 1 < — — (s) for < s < t. Now a change 
of variable in p will complete the theorem: 



p(t, x) 



/(0,X(0),y(0)) dv 
f (X(0),V(0))dVi 



1^(0) 



dvi 



i(0) cfo 



< 



0Vi 



F (Vi(0 J t,a,«),«a)- 5 -^(0) duidu 2 
Fq(w) dw. 



□ 



Proof of Theorems 1.3 and 1.5. We will first prove the result for (VP) and then 
comment on the necessary changes to adapt the proof to (RVP). Since we have 
defined the characteristics for (RVP) in (1), let us do the same for (VP), noting 

that the only difference occurs in the equation for — — : 

ds 

dV 

-j± = E 1 { S ,X{s)) 
V 2 {s,t,x,v)=v 2 > 

X(t, t, x, v) = x 
Vi(t,t,x,v) = vi 

Now, by the definition of E\ in (VP) and the well-known conservation of charge 
identity we find \E x {t, x)\ < \M for every ieK,i>0 and lim EUt, x) = ±-M. 
Using (12) and the compact support of fo in v%, 



\Vi(t)\ 



Vr(0) 



£i(s,X(s)) dx 



<C 



1 



-Alt 



for any characteristics (X(t),V(t)) along which f(t,X(t),V(i)) ^ 0. Taking the 
suprcmum over all such characteristics, we find 

Qi(t) < Ct 

for t large enough. The result can be seen to hold for (RVP), as well, since the 

equation for — — in (1) does not differ from that of (12) and the same bound 
ds 

\Ei(t,x)\ < \M again follows from charge conservation. 



To obtain the contrasting inequality for (VP), multiply the Vlasov equation by 
V\E\ and integrate over v to find 



Jv 1 E 1 d t fdv + Jv 2 1 E 1 8 x fdv + J Vl Eld v J dv 



0. 
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which is equivalent to 

d t J v x E x f dv - {dtEx) J vif dv + d x J vlE^fdv-p J i$f dv = J Effdv. (13) 

Let ji(t, x) = J vif(t, x, v)dv. As usual, integrating the Vlasov equation of (VP) in 
v yields the continuity equation 

d t p + d x jx = 0. 

Using this, we see that dtE\ = —j\ and integrating (13) in x yields 
d 



dt 



J I ViExf dvdx = J J pv\f dvdx - J jf dx + J J Ef,_ 



f dvdx 



d 
dt 



v\Eif dvdx 



p I vff dv — j\ ) dx + pE 1 dx 



(14) 



Now, by Cauchy-Schwarz the first term on the right side of (14) is nonnegative since 



Vif dv 



< 



v x f dv 



1/2 



f dv 



1/2 



and thus 



p I vjdv 



f dv 



vff dv)>[ / vif dv 



Jl 



The second term on the right side of (14) can be simplified as well 

-Mf\ M 3 



J P E\ dx=\J d x (Ef) dx = U^- 
Thus, from (14) we deduce 



12 



M 3 

— / / V\fEi dvdx > — — 



dt 



12 



and by integrating in time and using charge conservation, i.e. ||/(t)||i = ||/o||i = M, 
M 3 



12 



-t < 



E\(t,x)v\f(t,x,v) dvdx 
1 



Ei(0, x)vifo(x, v) dvdx 



< C^l + -M 2 Q!(i) 

< C(l + QiW). 

Hence, for t large enough Q\(t) > Ct. For (RVP) the analogous claim follows 
immediately from (4) since 

Qi(t) > Px(t) > Pl (t) > Pl (0) + \-Mt. 



Hence, for t large this implies 



and the proof is complete. 



Qi(t) > Ct 



□ 
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3. Behavior of solutions to (RVPN). 

Proof of Theorem 1.4-- We first obtain a priori bounds along the light cone as in 
[7]. Put 

e(t,x)= f ^/l + \v\ 2 Y,U{t,x,v) dv + ^E 2 (t,x), 
m(t, x) = vi^2 x ' w ) dv - 

J a 

As is well known, the local conservation of energy can be obtained by summing all 
of the Vlasov equations and integrating in v to find 

d t e + d x m = 0. 

An integration over all x then yields a bound on the total energy 



e(t, x) dx = J e(0, x) dx < C. 

Here, we must point out that this identity, and those which follow, hold only under 
the assumption of neutrality (N). If the plasma is not neutral (as is the case for the 
(RVP) system), then E\(t,x) — ► ±§M as x — > ±oo for every t > 0, where M^O 
is the total charge. In these cases, E(t,-) £ L 2 (M.) and the total potential energy 
is infinite. Now, for T > wc integrate the local energy identity over the forward 
cone originating from the point (0, x) 



= 



x+t 



(dte + d y m) dydt 



10 Jx-t 

and use Green's Theorem to obtain 

T 



(e - m)(t, x + t) dt+ (e + m)(t, x - t) dt 



x+T 



e(T,y)dy. (15) 



x-T 



Since e±m = J (y/l + \v\ 2 ± «i) £ Q f a dv + ±Ef > and 

e(T, y) dy < J e(T, y) dy = J e(0, y) dy, 

we have a bound on each term on the left side of (15) for all T > 0. Hence, the 
integral of e ± m is bounded along the rays of the cone: 



and 



J o °° J (VT+W-vi)Y,f°< dv + \ E i 



(t, x + t) dt<C 



(t, x-t) dt <C 



(16) 



for any x £ ffiL Next, we use a lemma which allows us to utilize the cone estimate 
in a way that controls the u-integral of ^ f a . 



Lemma 3.1. Let 



<J-(t, x)= / (Vl + H 2 - «i) ^2f a (t,x,v) dv, 
o-+(t,x) = J (Vl + \v\ 2 + v^j y^/ a (t, x,v) dv, 
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Then, 



k(t,x) = [ ^l + \v\*J2f a (t,x,v) dv. 

J a 
J a 

j ^2 f adv ~ 3 V CT + fc - 



Proof. We will show the former inequality, as the proof of the latter is similar. Let 
R > and split the w-integral: 



fjZfadv = [ Y j f a dv+ [ V/ Q 

J - J\v\<R ~ J\v\>R ~ 

< / > fa dv + ; 



dv 



Consider \v\ < R and notice 



VI + \v\ 2 - VI 



l + vl 



> 



> 



> 



y/l + \v\ 2 +V! 
4 



1 + vl 



l 



1 



2VT+W 



So, we find 



Jj2f a dv < J f° ■ (Vl + W - y i) ( Vl + R 2 ) dv + 



Vi + R 2 



< 2y/l + R 2 ■ a_(t,x) + 



k(t, x) 

7TTW' 



If < cr_ < k take R = J 1 so that \j\ + R 2 = 



Ef a dv < 2 J — ■ (T— 
\ rr 



and 



3 V c_fc. 



Now, if k < cr_ then 

f^f a dv< f v'l + H 2 ^ dv = k < < 3V^fc- 

Finally, if <t_ = then since f a > for all a = 1 , . . . , N 

Jj2f« dv = ° = 3V^fc 



and the proof is complete. 



□ 
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Since we are now considering the relativistic, neutral system (RVPN), recall that 
the characteristic equations are similar to that of (1) but must depend upon the 
charge of the corresponding species: 



dV lc 



dX a 
ds 



V la (s) 



e QJ Ei(s,X(s)) 



(17) 



ds 

V 2a (s,t,x,v) = v 2 
X a (t, t, x,v) = x 

Vl a (t,t,X,v) = V\ 



Now, let a be given and choose characteristics (X(t), V(t)) as in (17) such that 
f a (t, X(t),V(t)) 0. For brevity, we suppress the characteristic dependence upon 
the choice of a. Let 

Co = sup{|ui| : 3t e [0, l],x,V2 G K with f a (t,x,v) ^ 0}. 

and suppose that t > and V\{t) > 2Cq. Define 

A = sup |r e (0,t] : Vi(a) > ^(t) Vs E [t - r,t] 



and notice that 



so that t — A > 1 and 



Define 



Then, for t - A < s < t, 



-(X(s)-X c (s)) 



V l (t-A)>-V 1 (t)>C 



V 1 (t-A) = ±V 1 (t). 



X c (s) = X(t) + s-t. 



l + V 2 (s) 2 



< 4 



< 



i + y 2 (o) 2 



c 



Hence, integrating in time over [t — A, t] we find 



\X(s)-X c (s)\< 



CA 



(18) 



(19) 
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Integrating (17) and using the cone estimate (16), we find 



Vi(t) = V 1 (t-A) + e a f EkfaXcia)) ds 

Jt-A 



t-A 

(E 1 (s,X(s))-E 1 (s,X c (s))) ds 

t-A 

1/2 



< Vi(t-A)+ max|e Q 

V a 
I t-A JX c (s) 



< VAt- 



A 1/2 (f E%(a,X e (s)) ds 

ft fX(s) f 

/ / 22\e a \f a (s,x,v) dvdxds 

Jt-A JXc{s) J „ 

ft fX(s) p 

(t - A) + CA 1 / 2 + C / J2f a (s,x,v) dvdxds 

Jt-A JXr(s) J Z 



I t-A JX c (s) 

We use Lemma 3.1 and (19) on the last term to find 

ft fX(s) f ft pX^+CAV^t)- 2 

I \ l!>2f a dvdxds < 3 / / \/o^k dxds. (20) 

Now, the well-known conservation of energy identity gives a uniform bound on 
J k(t, x) dx and thus 

t ,X c (s)+CAVi(()" 2 ft j. 

/ k(s,x) dxds < / / k(s,x) dxds 

i-AH c (s) Jt-A J 

< f C ds 

Jt-A 

= CA. (21) 
Using (16), the integral of ct_ can be estimated as 

rt f X c (s)+CAV 1 (t)~ 2 

er_ (s, x) dxds 

t-A JX c (s) 



ft fX(t)-t+cAVi(ty 

/ / <r_(s,y + s) dyds 

Jt-A J X(i)-t 

x(t)-t+cAVi(ty 2 f t 

er_ (s, y + s) dsdy 



ix{t)-t Jt-A 

fX^-t+CAV^t)" 2 

< Cdy 
Jx(t)-t 

= CAVi{t)- 2 . 

Combining this with (21), the use of Cauchy-Schwarz in (20) gives us 

' [jlfa dvdxdt < 3(CA) 1 / 2 {CAVi{t)- 2 ^ 1/2 ' r ' A 

I t-A Jx c {s) 

and finally, we have 

A 



f r S [J2f- dvdxdt ^ ^ca) 1 ' 2 (cm®-*) 1 ' 2 < c*-A- 

Jt-A Jxjs) J „ 

Vi(t) < Vi(t - A) + CA 1/2 + C- * 



vi(ty 
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Using (18) this becomes 

1 



2 Vi(t)<C 



Hence. 



and 



Finally, 



Vi{t). 

CV x {t) 2 < A 1/2 V 1 (t) + A 
C+^V,{t) 2 <(±V,{t) + A^ 



CVi{t)=\^Jc + \-\\v 1 {t)<& 1 ' 2 . 
By definition, A < t and thus 

Vxit) < Ct 1 ' 2 

as long as V\(t) > 2Co- Similar estimates may be derived if we consider V\{t) < 
—2 Co- Combining all cases, we find 

< max{2C ,Ct 1/2 } 

whence 

|Vi(t)| < 2C a + Ct 1 / 2 . 

Since a is arbitrary, this estimate holds for every /„ and taking the suprcmum over 
all such characteristics, we finally have 



Qi(t) < cVT+t 

completing the proof. □ 



4. Behavior of solutions to (VPN). 



Proof of Theorem 1.6. Since we are now considering the classical, neutral system 
(VPN), recall that the characteristic equations are similar to that of (12), but as 
for (RVPN) they now depend upon the charge of the corresponding species: 



dX a 
ds 



V la (s) 



(22) 



^ = e a E 1 (s,X(s)) 
as 

V 2a (s,t,x,v) = v 2 
X a (t, t, x,v) = x 

Vl a (t,t,X,v) = Vi 

Notice, the only difference between (22) and (17) is the absence of the "hat" in 

the equation for — - — . We begin by deriving a well-known bound which limits 
ds 
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the growth of spatial moments for the Vlasov-Poisson system. Using the Vlasov 
equation, we find 



^2 // x2 ^2 x > w ) dvdx = 2 fc II XVl ? a dvdx 

a a 

= 2 J J v\ f a dvdx 

a 

+ 2 J J x e a f a (t, x, v)E\ (t, x) dvdx. 

The first term on the right side of the equation is dominated by the kinetic energy 
fC(t) =11 \v\ 2 fa(t, x, v) dvdx. 

J J a 

The second term satisfies 

J J x e a f a (t, x, v)Ei (t, x) dvdx = J pxE\ dx 

= j xd x {^ E l S j dx 
= -\J Ef{t,x) dx. 

So, when the two terms are added, the result is uniformly bounded by energy 
conservation. Hence, we find 



^2 j j x 2 ^2fa(t,x,v) dvdx 



< C 



and integrating twice in time 

^^2f a (t,x,v) dvdx < C(l + t 2 ). 

a. 

Using the form of the field, we find for x > 



\Ei(t,x)\ 



P{t, y) dy 



< c(max|e Q |) J J ^ f a (t,y,v) dvdy 

r°° r ^2 
- C J J z2fa(t-,y,v)-^ dvdy 



C 



1 + i 2 



This inequality is similarly obtained for x < so that 



\Ei{t,x)\ < Cmin 1, 



Now, using a method of Horst [11] we estimate characteristics. Let a be given 
and (X(t),V(t)) be characteristics such that f a (t, X(t),V(t)) ^ 0. Suppose that 
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t > and V\{t) > 0, and define 

t* = inf{r e [CM) : Vi(s) > for s G [r,i]}. 
Then, either Vi(r*) = or r* = and in both cases, V?{t*) < C. Using (22) 



V?(t) < F 1 2 (r*) + 2(ma X |e Q |)y |£i(s,X(s))| Fr(s) 

< c + c J r min ( lj ~^y^) ^ (s) ds 

< C + C J ndn(l,^^jx(s) ds 



ds 



Hence, 



X(t) 



< C + C I min ( 1,^- I dx 

Jx(t-) x 2 



C + C I min [ l,^-t_L I rf x 



'1+t 2 i-oo i i .2 > 

< C + C[ I 1 dx + I _ dx 

/i+t 2 x ; 



< c(l + Vl + i 2 ) 

< C(l + i) 



vi(t) < cvm. 

Similar estimates hold for V\{t) < and for characteristics of any f a since a was 
arbitrary. With this, we take the supremum over all such characteristics and obtain 
Qi(t) < C*v / TTI □ 
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